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Starting with the basic control system model often employed in NMR pulse design, we derive 
more realistic control system models taking into account effects such as off-resonant excitation 
for systems with fixed inter-qubit coupling controlled by globally applied electromagnetic 
fields, as well as for systems controlled by a combination of a global fields and local control 
electrodes. For both models optimal control is used to find controls that implement a set of 
two- and three-qubit gates with fidelity > 99.99%. 
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1. Introduction 

Controlling the dynamics of quantum systems is a problem of great recent interest with 
many applications from quantum chemistry to quantum computing. Since the dynamics of a 
quantum system is essentially determined by its Hamiltonian, a central theme in quantum 
control is Hamiltonian engineering. This is typically accomplished by the application of 
external fields, whose interaction with the system modifies its effective Hamiltonian, thereby 
changing its dynamical evolution. The type of control fields that are available tends to vary 
significantly depending on the type of system and the actuators available. In traditional 
applications such as nuclear spin engineering and nuclear magnetic resonance (NMR), the 
fields are often radiofrequency pulses, while in applications involving transitions between 
electronic states in atoms or vibrational modes of molecules, the control fields arc usually 
laser pulses. In more recent applications such as various solid-state architectures based on 
semi-conductor quantum dots or Joscphson junctions, the controls arc often electric fields 
produced by external electrodes. Despite the disparate physical realizations, all of these 
applications share important similarities, in that the external fields modify the Hamiltonian 
in a certain way, and the objective generally is to find a particular control that modifies the 
dynamics in a certain way, and achieves a certain objective, e.g., of steering the system to a 
particular state or realizing a desired quantum process. 

There are many differences in the specifics, however, which can be very important for the 
design of effective and experimentally feasible controls. In NMR applications, for example, 
we can often approximate the system well by assuming that the control fields induce rapid 
local rotations of individual spins, which interact through a fixed, weak coupling between 
them that is not affected by the control fields. Furthermore, we can often assume that we 
can selectively address individual spins by frequency-selective control pulses, and that we 
can perform local rotations about orthogonal axes. The same approximations, however, are 
not necessarily valid for other systems such as solid-state architectures. In particular for 
systems controlled by local voltage gates, it is often difficult to perform rotations about 
orthogonal axes, and a voltage applied to a local electrode may have non-negligible effects 
on the dynamics of nearby quantum dots, and in the presence of multiple control voltages 
there can be signficant cross-talk effects, all of which can signficantly complicate effective 
control design. In this paper we attempt to incorporate some of these effects to derive more 
realistic control system models for two types of systems of practical interest, and apply 
optimal control theory to find control solutions to implement a set of elementary gates for 
the resulting more complex control system models. 
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Basic NMR control system model and geometric control 

The starting point for many optimal control problems in quantum computing such as 
optimal implementation of quantum gates has often been a control system model inspired 
by liquid-state nuclear magnetic resonance (NMR). Using the common shorthand o'^"'' to 
denote an Af-factor tensor product in which the nth factor is the Pauli matrix aj^ for 
k S {x, y, 2} and all other factors are the identity I, the Hamiltonian for a system consisting 
of A'^ spin- 1 nuclei (qubits) in this model is usually given by 

JV 

n<n' n — 1 

where the first term is a fixed inter-qubit coupling term, in this case of Ising type, and we 
have 2N independent control fields U2n-l{t) and U2n{t) capable of performing local x and y 
rotations on the nth spin, respectively. 

If the interqubit couplings J„n' controllable, or at least can be switched off, then it is 
very easy to see that we can implement any local unitary operation U G SU{2) on the nth 

(n) 

qubit, e.g., by performing a sequence of three local rotations about the orthogonal axes iri 
and o'y"'', e.g., using the standard Euler decomposition U = Ux{ot)Uy{l3)Ux{l) with 
Ux{oi) = exp(— iao-j:), etc, for suitable values of the rotation angles a,f},^, and two-qubit 
gates can be implemented using the Cartan decomposition ifll [3l 

(7 = C/i [Uy Z{a3)U-^] [C/-1 Z{a2) U-A Z(a^) U2, 

where Z{a) = exp(— iao-^^'o-^^' ) corresponds to free evolution under the Ising-coupling 
Hamiltonian for a suitable time, Oi and IJ2 correspond to local operations on both qubits, 
which can be implemented using the Euler decomposition, and Ux = Ux^\^) ui^\^) 

correspond to local x-rotations by ^ on both qubits, and similarly for Uy. 

If the inter-qubit coupling terms J„n' fixed and not controllable then this gives rise to 
a drift term, which complicates the situation. If the coupling is weak and we can apply 
"hard pulses", e.g., if we can ensure 

V n' 



then the drift term is negligible (as far as the implementation of local gates is concerned), 
and the geometric control pulse sequences derived from the Euler and Cartan decomposition 
are often a good approximation. However, if the coupling is stronger or hard pulses are not 
available, neglecting the drift term will result in inaccurate gates. Moreover, in realistic 
applications, hard pulses may lead to other problems such as off-resonant excitation. 

Although we can in principle try to find geometric decompositions that take the drift 
term into account (5), geometric control pulse sequences derived from Lie group 
decompositions are usually not optimal even when the drift is negligible, and optimal 
control has been shown to be able to produce more effective or efficient control pulses with 
regard to various perfomance indices such as gate operation timesor average fidelity even for 
cases where geometric decompositions are valid, at least theoretically [See e.g. iffit]. 



Realistic control system models for optimal control 

The success of optimal control theory in finding more effective control pulse sequences for 
various model problems has fueled a surge of interest in this area. However, the control 
system models on which many of the optimal control calculations so far have been based are 
generally too simple and often unrealisitc for the systems they are supposed to model. The 
model described in the previous section, for example, is based on many approximations that 
may be questionable even in NMR settings, and certainly for other applications. For 
instance, in the ideal case the control fields U2n—i{t) and U2n(i) induce only an x-rotation 
and j/-rotation, respectively, on the nth qubit, but the desire to speed up gate 
implementation times may produce strong control fields likely to induce off-resonant 
excitations not accounted for the in model. 

Off-resonant excitation means that the fields are not only felt by the target qubit but also 
by other qubits. In the simplest case we can try to model off-resonant excitation by 

(n) (n) 

replacing the individual interaction Hamiltonians and a-y ' for the independent control 
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fields M2n-i{t) and U2n{t) by linear combinations of the coupling terms 

n' 

H2n = "2n,2n'-l'5'i"' 



' I ^ /-/ \ /-II _ _ . . /T ^ ' _1_ <-V „ „ , /T ^ ' 



The new interaction Hamiltonians Hm for m = 1, . . . , 2Af no longer correspond to local 
rotations about mutually orthogonal axes, which makes finding geometric decompositions 
for this case very challenging in general. However, even arbitrary (linear) cross-talk need not 
be an obstacle for optimal control as shown in Q for a model system of this type. Provided 
the cross-talk Hamiltonians are known, based on a combination of theoretical models and 
experimental characterization, it still appears to easy to find effective control pulses using 
optimal control algorithms. The main drawback is that this simple model does not 
accurately model real cross-talk. 



Control via globally applied fields with off-resonant excitation 

The intrisic system Hamiltonian for a system of N spins with resonance frequencies ui„ with 
static Ising coupling J„„i in a fixed laboratory frame is 

Ho = h Jnn'f'z - h2_^ —a\ ' . 

n<n' n 

If we apply a global electromagntic field B(t) with components Bx(t) = Bit) cos(t<jt + cji) and 
By{t) = Bit) sin((i;t -f 0), respectively, the resulting control Hamiltonian is 

He = - | ^7uS(t)[cos(a.;t + - sin(Ljt + 0)4"'] 



Transforming to a multiply rotating frame defined by Us(i) = exp{— j| X]n '^"''^z"')i the 
interaction picture rotating wave approximation Hamiltonian is {7) 

H = hJ2 -4"'' - ^^7n-B(t)[cos(Aa;„t + 0)4"' -sin(Aa;„t + <^)4"'] 



where Auin =10 — LOn and f2n(i) = 7nB(f). If the field is exactly resonant with the nth spin, 
uj = uJn, then the term corresponding to the nth spin simplifies to 

-^7nB(t)[cos(0)4"' - sin(0)4"']. 

If we can further argue that all other spins have frequencies Wn far detuned from ui, \Bit)\ is 
slowly varying and not too large, |7„/B(t)| <^ |Aijj„/| for all n' n, then for sufficiently 
large T the rapidly oscillating terms will "average to zero", i.e., we can assume 



1 

- / Bit)-f„,cosiAui„,t + <t,)dt^O 
Jo 

1 

- B(t)7„'sin(Aa;„,t + 0)a!tRiO. 

In this case the effective control Hamiltonian associated with the field B(t) simplifies to 



= -^n„(t)[cos(</.)4"' - sin(0)4"']- 



Thus, if we are able to simultaneously apply TV independent fields i?i"\t) cos(aJnt) with 

frequencies uJn and A'^ independent fields By"\t) siniuint) with frequencies oj„, for instani 
then assuming all the previous assumptions and approximations still hold, we obtain 
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exactly the Hamiltonian with M2n-i (i) = -|7nSi"\t) and 

In practice, it is questionable, however, whether these approximations are still valid in 
general, especially if we simultaneously apply 2N control fields _Bi"'{t) and By"'\t), all of 
which may be strongly modulated in time as a result of an optimization procedure, and thus 
not slowly varying. Furthermore, the interaction picture RWA Hamiltonian |6]l is not a very 
good starting point for a general optimization either because each of the independent 
control Hamiltonians 

H,n(t) = - ^ 7n [cos( Aa;„,^f + (l>,n)ai"^ - sin(Atj„,mi + </>m)CT^"'] 
n 

with LOn^rn = — LOn for a field Bm(t) with Bx{i) = Bm{t) cos{LUmt + (pm) and 
By{t) = Bm(t) sin{uimt + 0m)i in the total Hamiltonian 

m 

is time-dependent, containing rapidly oscillating sine and cosine terms. Therefore, wc might 
as well let the component functions Bx{t) and By(t) of the field be completely arbitrary and 
work with the general Hamiltonian 

n n 

with Hq as in Eq. (Ull, instead. 



Control via fixed global fields and local control electrodes 

In some applications such as NMR the only type of control we have is globally applied 
control fields. However, for many systems, especially solid-state systems such as the nuclear 
spins of donor atoms embedded in a substrate , or electron spins in quantum dots Ql , for 
example, some local control in the form of control electrodes is available, and in a typical 
setting, control is achieved by combining globally applied electromagnetic fields with local 
fields provides by control electrodes (| jQj I. In this case we often have a global electromagnetic 
field B(t) with components Bx{t) = B{t) cos(a;4), By(t) = B{t) sm{uit), in addition to a 
static magnetic field in the z direction, and a Hamiltonian of the form 

n<n' " 

- ^ J]]7nB(t)[cos(aJt)o-i"' ~ sm{LL)t)a^y"^], 

n 

as before, but in addition we now have some degree of control over the resonance frequency 
LU„ of the nth qubit via local control electrodes. Note that unlike in the previous model we 
have assumed isotropic Heisenberg coupling, i.e., 

• = ai"^a^ ' -|- d-y"^ Uy" ^ + ai"^ ai" \ instead of ai"^ai" •' Ising coupling as this 

seems to be more appropriate for solid-state systems, although the assumption of isotropy 
may not be appropriate for some systems. 

One control approach is to simultaneously vary the voltages applied to the local control 
electrodes as well as the globally applied field. However, it practice it is much simpler to 
apply a global field B{t) with a fixed amplitude B and fixed frequency u! and only vary the 
voltages applied to the local control electrodes. In this case we can introduce a rotating 

frame, this time rotating at the field frequency uj, Us{t) = exp(— i| Xln'^^z"')' *° obtain 
the slightly simpler interaction picture Hamiltonian 

7i<n' n n 

where we set SI = |B7o and take 7„ to be the coupling constants 7„ in units of some 
coupling strength 70, so that for a homogeneous system consisting of a chain of spins that 
couple identically to the field we have 7,1 = 1 for n > 0. 

In this model the constants J„„' and 7„ are usually fixed and determined by the physics 
of the system. The frequency uj and Rabi frequency Q of the field can be chosen initially but 
are then fixed. Thus, all the dynamic control is achieved by dynamically changing the 
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control voltages and thus a;„{V). This type of control therefore leads to a different 
partitioning of the total Hamiltonian into a fixed drift and variable control part 



In general, due to cross-talk effects flTJ the energy level shift of the nth qubit may be a 
complicated function of all the voltages applied, not only to the electrode directly above it, 
but to neighbouring control electrodes as well, i.e., LOn may be a function of all the control 
voltages V = (Vi, . . . , Vm) applied, i.e., uJri = Wn(V). 



4. Optimal gate implementation 
4.1. Optimal control formulation 

One way to solve the problem of finding controls u(t) that implement a desired set of 
quantum gates, especially for control systems that are not amenable to simple geometric 
decomposition schemes, such as the ones discussed in the previous section, is to formulate 
the task as an optimization problem. This typically involves (a) choosing a target function 
to be optimized, (b) finding a suitable parametrization for the set of admissible (and 
experimentally feasible) controls, and (c) finding a solution to the resulting, possibly 
constrained, optimization problem based on application of Pontryagin's maximum 
principle il2t . 

For the task of implementation of quantum gates or processes, a natural choice for the 
objective function is the gate fidelity l| jgjl 

r{UT,U) = ^'RTv{U^U). (13) 

The reason for the factor of l/N, where A'^ is the system dimension, is to ensure that 
varies between —1 (for U = —Ot) and +1 (for U = Ut) independent of the system 
dimension. In general, maximizing the gate fidelity is equivalent to minmizing the gate error 

e(UT,u) = \\u -ijTf. (14) 

This is very to see if the target operator Ut and the actual operator U are both unitary 
(and thus (/''U = IJIUt = /) as 

lit/ - UtW^ = Tt[{U - Ut)Hu - Ut)] 

= tt[u^u + u\,Ut - (i^iiT - ulu\ 

= 2Ar- 2»Tr(l/^(7). 

Noting that the implemented operator U = U,i{tp) is a function of the control u, the control 
objective is to find an admissible control Uopt such that 

^{Ut, U^apt (tp)) = ma^HUT, U^itp)), (15) 

ug.4 

and of course, Uu{t) has to satisfy the Schrodinger equation 

^Mt) = ^lHHt)]U^{t), (16) 
dt n 

where H[u{t)] is the control-dependent Hamiltonian of the system. 



4.2. Iterative control optimization algorithm 



Once the Hamiltonian H[u(t)], the target time tp and the target operator Ut have been 
chosen, the fidelity depends only on the choice of the control u(t), i.e., = jFu Now let u 
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and u + Au be two controls. Noting tfiat 

J to 

and that for a control-linear Hamiltonian H[u{t)] = Hq + Um(t)Hrn we have 



)dT 



we obtain immediately 



-5RTr[C/T{(7„+Au(i)-C/u(t))] 



Nh 
2 



-STr 



if" U^{tF,t)J2Au^{t)H^U^+A,,{t)dt 
J2 [ " Aurn{t)^Tr\ulu^(tF,t)HmU^+AM] dt 



which shows that setting 



Au,„{t) = e„(t)STr |C/t ;7u{tF,t)Jfmi7u+Au(t) 



with €m{t) > for all t will increase the fidelity. Thus, starting with any initial guess uo(t), 
iteratively solving the Schrodinger equation while updating the control in e ach iteration 
according to the rule u("+^) = u'"' + Au'"', with Au'"' chosen as in Eq. | |l9| , will 
monotonically increase the fidelity in each iteration, and as the fidelity is of course bounded, 
the iteration must converge to some limiting fidelity J^°° , although we cannot guarantee 
that the limiting value will be equal to the global maximum of the fidelity. In fact, many 
optimal control algorithms are based precisely on this iterative update scheme. 



Parametrization of controls and discretization 

To solve the Schrodinger equation numerically, we still require some discretization. The 
simplest approach is to assume piecewise constant fields u(t) = \if^ for t^—i <t< t^, in 
which case the exact solution of the Schrodinger equation is explicitly 

U^tp) = {/u^(At^)i7u^_i(AtK_i)...l/ui(Aii) 

where Atf; = tf^ — tfc_i, tx = tp, and 



-AtkH{uk) 
n 



U^^(Atk) = exp 

is a simple matrix exponential and the update rule becomes 

r'rt fl 



Aum(tk) = em(t)S Tr [t/j,(7u^ (tx) ■ ■ ■ t^u^+i (tfc+i)//„C/u JAtfe) ■ ■ ■ (7ui {Ati)J . 

With this type of discretization there are two ways of implementing the update rule. In the 
global update approach all the matrix exponentials for all times tj. are evaluated to solve 
the Schrodinger equation, and the fields are updated for all times tk simultaneously at the 
end of each iteration step. This is essentially the well-known GRAPE algorithm as described 
in ^14h An alternative is to update the fields Um{t]^) at each time step tj,. Our analysis 
above shows that both approaches will monotonically increase the fidelity provided 
^m(tk) > for all t]^, although they will tend to converge to different solutions for the fields. 

In our limited numerical experiments for the systems considered here, both approaches 
seemed equally effective in finding satisfactory control fields but the local update approach 
tended to converge faster in terms of the computational time required to find a solution. 
The main reason for this appears to be that in the local update approach only a single Ufc is 
incremented at a time, and hence we can optimize the multiplier with very little 
computational overhead to speed up the convergence, as any local change of requires only 



Journal of Modern Optics 



2008JMO 



re-evaluation of a single matrix exponential (and two matrix multiplications and a trace 
operation) to assess its effect on the fidelity, while in the global update any change in e 
requires the re-evaluation of all the K matrix exponentials and K — 1 matrix multiplications 
to assess the fidelity of the new control. An alternative to this simple step-size control is to 
use higher derivatives to predict the optimal step-size but this also involves potentially 
significant computational overhead to find the best step size, in addition to the need to 
calculate analytical expressions for these derivates, thus it is not clear if this will yield a 
computational advantage or not, although it is worth investigating. 

A general shortcoming of this simple optimization procedure is that it does not take 
constraints into account or include penalities for less desirable fields. Simple constraints such 
as upper bounds on the field amplitudes can easily be incorporated into the algorithm. More 
complex constraints such as limiting the frequency bandwidth of pulses or minimizing the 
overall pulse energies etc require augmentation of the objective functional by constraint or 
penalty terms and modifications to the update rules. For some types of penalty terms 
iterative algorithms very similar to the algorithm above exist [see e.g. 03] i but a detailed 
discussion of appropriate penalty functions, resulting iterative update rules, and comparison 
of the resulting optimal fields is beyond the scope of this paper. 



Simulation details and results 

The specific task we wish to accomplish here is finding controls u(i) that implement the 
following complete set of universal gates (i.e., target operators) 

C/t 6 {I I, Had (g) I, T I, I Had, I ® T, CNOT} 

for a two-qubit system, where we employ the common appreviations Had = exp(i^y), 
T = exp(i|Z) and CNOT = e"*'"''"' diag(7', X). The phase factor for the CNOT is necessary 
to ensure that the target operators are in (3it(4) since our model Hamiltonians have zero 
trace and hence can only generate operators in SH(4). As a second, more challenging task 
we consider the implementation of the Toffoli-like three-qubit gate 

IJt = diag(l, 1, 1, 1, 1, 1, iX) e Sil(8). 

For the simulations we assume uniform nearest-neighbour coupling, i.e., J„„/ = J&n.n' — i- 
Since the coupling constants J in our model have units of frequency, by choosing time in 
units of 1/J, we can set the coupling constants to unity. All other frequencies, unless 
otherwise states, are also chosen in units of J. We can furthermore choose the units of 
energy such that h, which has units of energy X time, equals 1, and can thus be omitted. For 
the two-qubit gates we set the target gate operation time to tp = 1, while for the three-qubit 
gate we choose tj? = 5 as a Toffoli gate normally requires five two-qubit gates. We choose 
these values as reasonable target times although we were able to find solutions for smaller 
tp for many gates, suggesting that the gate operation times could in fact be improved. 



Global control fields with off-resonant excitation 

For our first model system with Hamiltonian ( |10| , assuming 7„ = 70 for all n, and setting 
u\{i) = — |7o-Ba;(t) and U2{t) = |7o-By(t), the Hamiltonian becomes 

H = ZZ - i(a;iZI -(-a;2lZ) + Mi(t)(XI + IX) + M2(t)(YI + lY) 

for a two-qubit system, and for a three-qubit system we have 

H =(ZZI -I- IZZ) - i(a;iZn + W2IZI -I- W3IIZ) 

+ Mi(t)(Xn + IXI + IIX) + ui{t){Yll + lYI + IIY). 

Fig.^shows the actual control field solutions for a system with ui-^ = 10 and L02 = 12, i.e., 
a system with Ising coupling frequency about 10% of the smaller qubit frequency, i.e., strong 
enough to be non-negligible for single qubit operations, and a qubit frequency difference 
large enough to allow some selectivity but still small enough so that off-resonant excitation 
is a concern for pulses of the magnitude and duration required for to achieve the desired 
gate operation times. For this model system we choose a relatively large number of time 
steps, K = 100, motivated by the assumption that the x and y components of the globally 
applied electromagnetic field could be varied more or less arbitarily, although simulations 
suggest that we can find piecewise constant solutions with much fewer steps, if necessary. 
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0.5 

t (units of 1/J) 



0.5 

t (units of 1/J) 



0.5 1 

t (units of 1/J) 



Figure 1. Control field solutions for system | |25| l witii = 10 and 0J2 = 12 for elementary gates | |23| l and 
tf = 1 witii 100 time steps and no field constraints. 

The solutions are not unique and in fact different initial fields or update rules generally 
produce different solutions. Again, despite the increased complexity of the model 
(off-resonant excitation, no RWA) and the reduced number of control fields (two vs four in 
the simpler model Q), it appeared easy to find solutions that achieved the target fidelity of 
> 0.9999, or gate errors < 0.01%. However, as the figure shows, the magnitudes of the fields 
for at least some of the gates are relatively large. Further analysis also shows that the 
spectral range of the optimal controls produced by this algorithm tends to be rather broad, 
which may be a concern in practice when there are bandwidth limitations as filtering of high 
or low frequency components will reduce the gate fidelity. Some numerical simulations for 
related problems suggest that solutions with much lower spectral bandwidth often exist but 
it appears that a different parametrization of the fields and optimization approach are 
required to find such solutions. Finding controls to implement the three-qubit gate proved 
more of a challenge for this system, in part because convergence was extremely slow. 
Fig. ??(left) shows a solution with fidelity > 99.99%. Notice that the fields have high 
amplitudes and are very noisy. Fig. |4rieft) also shows that the fields have very large spectral 
bandwidth as expected, which woula make faithful implementation in the laboratory 
challenging. 



5.2. Control via fixed global fields and local control electrodes 

For our second model system with Hamiltonian \12\ , assuming again 7„ = 70 for all n, and 
setting Un(t) = — |aj„(V(t)) for n = 1, . . . , N , the Hamiltonian becomes 



H = XX + YY + ZZ - n{Xl + IX) + ui{t)Zl + U2it)IZ 
for a two-qubit system, and for a three-qubit system we have 



(27) 



H =(XXI + YYI -I- ZZI -I- IXX -I- lYY + IZZ) - n{XU + IXI + IIX) 
+ ui(t)ZU + M2(t)IZI + M3IIZ. 



(28) 



Notice that we take u„{t) to be independent controls here to circumvent the complex and 
highly architecture specific problem of mapping the actual control voltages V onto the 
frequency shifts lu„(V), although it must be stressed that optimal control solutions for u(i) 
will only be useful for a specific system if the u(t) can be actually realized through proper 
choice of the voltages, which is likely to impose additional, non-trivial constraints. 

Again, we chose target gate operation times of = 1 and tp = 5 for the two- and 
three-qubit gates, respectively, but in this case we tried to minimize the number of steps K, 
assuming that frequent switching of the control voltages would be experimentally more 
challenging and increase errors. As Fig. [2] shows for f! = 10 (in units of J) we were able to 
find satisfactory solutions for all six two-qubit gates with only K = 10 switches. We note 
that f2 Si 10 appeared to be optimal in the sense that we had to increase K significantly for 
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U^=l ® I, fid=0.9999 



U^=T®l,fid=0.999901 
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Figure 2. Control field solutions for system |27| with £7 = 10 for elementary gates l |23[ l and tf = 1 with 
10 time steps and no field constraints. 
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Figure 3. Control field solutions for three-qubit Toffoli-like gate for model system ( |26[ l with 
UJ2 = 12 and 0^3 = 8 (left), and model system l|28| with f! = 10 (right). 



= 10, 



larger or smaller values of Q. For f! = 1, for example, we had to increase the number of steps 
to if = 40 and double the target time to be able to implement the first five two-qubit gates 
(corresponding to local operations) with fidelity > 0.9999, and for the final CNOT gate the 
fidelity we achieved was still only 0.9996, i.e., below threshold. For £7 = 5 we were able to 
find solutions above threshold fidelity for all target gates for K = 20; for K = 10 the fidelity 
of the CNOT gate remained below threshold (99.33%). Similarly, for f2 = 20 we were able to 
find solutions above the threshold fidelity for all gates for K = 20 but not for K = 10. For 
even larger Q, such as 30 or 40, K had to be increased to 30, and for Q in the 50-60 range 
K = 50 steps were required. These results can be partially explained considering that if f2 is 
on the order of the fixed coupling strengths then the weight of the drift term makes it 
harder to implement local operations hence the increase in the time required. For large fl 
local operations become easier to implement but entangling gates such as the CNOT, which 
rely on generation of entanglment via the fixed Ising coupling term, become more 
challenging, requiring more complex pulse shapes, hence necessitating a larger number of 
steps if. It is interesting to note that for this system we had no trouble finding solutions to 
implement the three-qubit gate with fidelity > 99.99%. A possible solution is shown in 
Fig. Isiright). Although the pulse profiles appear more complicated, their magnitudes are 
actually lower than for some of the two-qubit gates, and although K = 50, the number of 
switches per time unit is still only 10. Unlike for the previous model Fig. WIright) shows that 
the spectral profile and bandwidth of the solutions are much more desirable. 



May 19, 2008 



16:44 Journal of Modern Optics 



2008JMO 



10 



REFERENCES 




Spectral decomposition of fields (b) 
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Figure 4. Spectrum of control fields for three-qubit Toffoli-like gate for model system l |26| with uj\ = 10, 
u)2 = 12 and 1^3 = 8 (left), and model system \28\ with f! = 10 (right). 

6. Conclusion 

Starting with the most basic control system model commonly used in applications such as 
(liquid-state) NMR, we have derived more realistic control system models for two types of 
control of special interest for many applications: (a) control of qubits with fixed inter-qubit 
coupling using globally applied electromagnetic fields in the regime where off-resonant 
excitation by the control fields is non-negligible, and (b) control of qubits with fixed 
coupling using a combination of a fixed globally applied control field and local control 
electrodes to shift the resonance frequencies of the qubits. Despite the substantially 
increased complexity of the resulting models, reduced number of independent controls, and 
non-trivial constraints on the allowed frequency of voltage changes in the latter model, 
optimal control techniques even in the rather basic form implemented, overall still proved 
very effective in finding solutions that achieved the target fidelity of 99.99% for almost all 
cases, although there appears to be scope for further improving the controls with regard to 
amplitude and bandwidth considerations, and for improving the overall convergence of the 
algorithm, which even with some of the modifications discussed, appears to be rather slow 
especially for the three-qubit gates. A somewhat unexpected conclusion of the calculations is 
that our second model system, which is controlled mainly by local voltage gates that can 
only alter the detuning of the qubits from the fixed (unmodulated) global field appears to be 
easier to control, even when constraints on the number of voltage changes are imposed, than 
the first model system, which is controlled by global fields with arbitrary pulse envelopes. 
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